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ABSTRACT 



38 We present multiwavelength studies of the 106.6 ms 7-ray pulsar PSR 

39 J1907+06 near the TeV source MGRO J1908+06. Timing observations with 

40 Fermi result in a precise position determination for the pulsar of R.A. = 

41 19 h 07 m 54!7(2), decl. = +06°02'16(2)" placing the pulsar firmly within the TeV 

42 source extent, suggesting the TeV source is the pulsar wind nebula of PSR 

43 J 1907+0602. Pulsed 7-ray emission is clearly visible at energies from 100 MeV to 

44 above 10 GeV. The phase-averaged power-law index in the energy range E > 0.1 

45 GeV is T = 1.76 ± 0.05 with an exponential cutoff energy E c = 3.6 ± 0.5 GeV. 

46 We present the energy-dependent 7-ray pulsed light curve as well as limits on off- 

47 pulse emission associated with the TeV source. We also report the detection of 
4B very faint (flux density of ~ 3.4 /xJy) radio pulsations with the Arecibo telescope 

49 at 1.5 GHz having a dispersion measure DM = 82.1 ± 1.1 cm~ 3 pc. This indi- 

50 cates a distance of 3.2 ± 0.6 kpc and a pseudo-luminosity of L1400 — 0.035 mJy 

51 kpc 2 . A Chandra ACIS observation revealed an absorbed, possibly extended, 

52 compact ( 4") X-ray source with significant non-thermal emission at R.A. = 

53 19 h 07 m 54!76, decl. = +06°02'14.6" with a flux of 2.3±?;| x lO-^ergcm-V 1 . 

54 From archival ASCA observations, we place upper limits on any arcminute scale 

55 2-10 keV X-ray emission of ~ 1 x 10 _13 erg cm _2 s _1 . The implied distance to the 

56 pulsar is compatible with that of the supernova remnant G40.5— 0.5, located on 

57 the far side of the TeV nebula from PSR J1907+0602, and the S74 molecular 
5B cloud on the nearer side which we discuss as potential birth sites. 

59 Subject headings: pulsars: individual: PSR J1907+0602 — gamma rays: obser- 

60 vat ions 



61 1. Introduction 

62 The TeV source MGRO J1908+06 was discovered by the Milagr q collaboration at a me- 

63 dian energy of 20 TeV in their survey of the northern Galactic Plane (lAbdo et al.ll2007l ) with 

64 a flux ~ 80% of the Crab a t these energies. It wa s subsequently d etected in t he 300 GeV - 

65 20 TeV range by the HESS f Aharonian et al.ll20Q9h and VERITAS JWardll2008h experiments, 
ee The HESS observations show the source HESS J1908+063 to be clearly extended, spanning 

^0.3° of a degree on the sky with hints of energy-dependent substructure. A decade earlier 



Lamb fc Macomb! ( 119971 ) cataloged a bright source of GeV emission from the EGRET data, 



69 GeV J1907+0557, which is positionally consistent with MGRO J1908+06. It is near, but 
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i nconsistent with, th e third EGRET catalog (IHartman et a 



. 1999) source 3EG J1903+0550 



([Roberts et al.l l200lh . The Large Area Telescope (LAT) flAtwood et al.l 120091 ) aboard the 
Fermi Gamma-Ray Space Telescope has been operating in survey mode since soon after its 
launch on 2008 Jun e 11, carrying out continuous observations of the GeV sky. The Fermi 
Bright Source List (jAbdo et al.ll2009bl ). based on 3 months of survey data, contains 0FGL 
J1907.5+0602 which is coincident with GeV J1907+0557. The 3EG J1903+0550 source loca- 
tion confidence contour stretches between 0FGL J1907. 5+0602 and the nearby source 0FGL 
J1900. 0+0356, suggesting it was a conflation of the two sources. 

The Fermi LAT collaboration recently reported the discovery of 16 previously-unknown 
pulsars by using a time differencing technique on the LAT photon data above 300 MeV 
JAbdo et al.l l2009at ). 0FGL J1907.5+0602 was found to pulse with a period of 106.6 ms, 
have a spin-down energy of ~ 2.8 x 10 36 erg s _1 , and was given a preliminary designation 
of PSR J1907+06. In this paper we derive a coherent timing solution using 14 months 
of data which yields a more precise position for the source, allowing detailed follow-up at 
other wavelengths, including the detection of radio pulsations using the Arecibo 305-m radio 
telescope. Energy resolved light curves, the pulsed spectrum, and off-pulse emission limits at 
the positions of the pulsar and PWN centroid are presented. We then report the detection of 
an X-ray counterpart with the Chandra X-ray Observatory and an upper limit from ASCA. 
Finally, we discuss the pulsar's relationship to the TeV source and to the potential birth 
sites SNR G40. 5-0.5 and the S74 Hll region. 



2. Gamma-ray Pulsar Timing and Localization 



91 Th e discovery and initial pulse timing of PSR J1907+06 was reported by lAbdo et al. 

92 (j2009al ). The source position used in that analysis (R.A. = 286.965°, Decl. = 6.022°) was 

93 derived from an analysis of the measured directions of LAT-detected photons in the on-pulse 

94 phase interval from observations made from 2008 August 4 through December 25. Here, we 

95 make use of a longer span of data and also apply improved analysis methods to derive an 

96 improved timing ephemeris for the pulsar as well as a more accurate source position. 

97 For the timing and localization analysis, we select ed "diffuse" class ph otons (events 

98 that passed the tightest background rejection criteria (lAtwood et al.l 120091 )) with zenith 

99 angle < 105° as is standard practice and chose the minimum energy and extraction radius 

100 to optimize the significance of pulsations. We accepted photons with E > 200 MeV from 

101 within a radius of 0.7° of the nominal source direction. We corrected these photon arrival 
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102 times to terrestrial time (TT) at the geocenter using the LAT Science Tool@ gtbary in its 

103 geocenter mode. 
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We fitted a timing model using TEMPo2 (jHobbs et al.ll2006l ) to 23 pulse times of arrival 
(TOAs) covering the interval 2008 June 30 to 2009 September 18. We note that during 
the on-orbit checkout period (before 2008 August 4) several instrument configurations were 
tested that affected the energy resolution and event reconstruction but had no effect on the 
LAT timing. To determine the TOAs, we generated pulse profiles by folding the photon times 
according to a provisional ephemeris using polynomial coefficients generated by Tempo2 in 
its predictive mode (assuming a fictitious observatory at the geocenter). The TOAs were 
measured by cross correlating each pulse pr ofile with a kern el density template that was 
derived from fitting the full mission dataset ( IRay et al.l 120091 ) . Finally, we fitted the TOAs 
to a timing model that included position, frequency, and frequency derivative. The resulting 
timing residuals are 0.4ms and are shown in Figured! The best-fit model is displayed in 
Table [TJ The numbers in parentheses are the errors in the last digit of the fitted parameters. 
The errors are statistical only, except for the position error, as described below. The derived 
parameters of E, B , and r c are essentially unchanged with respect to those reported by 



Abdo et al.l (j2009al ). but the position has moved by 1.2'. 



The statistical error on the position fit is < 1"; however, this is an underestimate of the 
true error. For example, with only one year of data, timing noise can perturb the position 
fit. We have performed a Monte Carlo analysis of these effects by simulating fake residuals 
using the FAKE plugin for Tempo2. We generated models with a range of frequency second 
derivatives (±2 x 10 _22 s~ 3 , the allowed magnitude for v in our fits) to simulate the effects 
of timing noise and fitted them to timing models. Based on these simulations, we assigned 
an additional systematic error on the position of 2", which we added in quadrature to the 
statistical error in Table [TJ As a result of the improved position estimate provided by this 
timing analysis, we have adopted a more precise name for the pulsar of PSR J1907+0602. 



3. Detection of Radio Pulsations 



129 To search for radio pulsations, we observed the timing position of PSR J1907+0602 

130 with the L-wide receiver on the Arecibo 305-m radio telescope. On 2009 August 21 we made 

131 a 55-minute pointing with center frequency 1.51 GHz and total bandwidth of 300 MHz, 

132 provided by three Wideband Arecibo Pulsar Processors (WAPPs, iDowd et al.l ( 120001 ) ) . each 

133 individually capable of processing 100 MHz. We divided this band into 512-channel spectra 



http://fermi.gsfc.nasa.gov/ssc/data/analysis/documentation/index.html 
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134 accumulated every 128 /xs. The small positional uncertainty of PSR J1907+0602 derived from 

135 the LAT timing means that a single Arecibo pointing covers the whole region of interest. 

136 After excising strong sources of radio-frequ ency interference wi th rf if ind, one of the 



137 routines of the PRESTO signal analysis package (IRansom et al.ll2002l ). we performed a search 

138 by folding the raw data with the Fermi timing model into 128-bin pulse profiles. We then 

139 used the PRESTO routine prepf old to search trial dispersion measures between and 1000 
wo pc cm~ 3 . We found a pulsed signal with a signal-to- noise ratio S/N = 9.$ and duty cycle 

141 of about 0.03 at a dispersion measure DM = 82.1 ± 1.1 cm -3 pc. This value was estimated 

142 by dividing the detection data into 3 sub-bands and making TOAs for each sub-band and 

143 fitting for the DM with TEMPO. 

144 We applied the same technique for 4 different time segments of 12.5 minutes each and 

145 created a time of arrival for each of them. We then estimated the barycentric periodicity of 

146 the detected signal from these times of arrival. This differs from the periodicity predicted 

147 by the LAT ephemeris for the time of the observation by (—0.000005 ± 0.000020) ms, i.e., 

148 the signals have the same periodicity. 

149 Subsequent radio observations showed that the phase of the radio pulses is exactly as 

150 predicted by the LAT ephemeris, apart from a constant phase offset (depicted in Figure [2]) 

151 A confirmation observation with twice the integration time (1.8 hr) was made on 2009 

152 September 4. The radio profile is shown in the bottom panel of Figure. [2] with an arbitrary 

153 intensity scale. The pulsar is again detected with S/N of 3.4, 5.1, 7.3 and 8.6 at 1170, 1410, 

154 1510 and 1610 MHz. The higher S/N at the higher frequenc ies suggest a posit ive spectral 



155 index, similar to what has been observed for PSR J1928+1746 (ICordes et al.ll2006l ). However, 

156 this might instead be due to scattering degrading the S/N at the lowest frequencies — for 

157 the band centered at 1610 MHz the pulse profile is distinctively narrower (about 2% at 50% 

158 power) than at 1410 or 1510 MHz (about 3%). At 1170 MHz the profile is barely detectable 

159 but very broad. This suggests an anomalously large scattering timescale for the DM of the 
wo pulsar. Observations at higher frequencies will settle the issue of the positive spectral index, 
lei For the 300 MHz centered at 1410 MHz, where the detection is clear, we obtain a total S/N 

162 of 12.4. 

163 With an antenna T sys = 33 K (given by the frequency-dependent antenna temperature of 



2 This was estimated using another software package, SIGPROC (a package developed by Duncan Lorimer, 
see http://sigproc.sourceforge.net/), which processes the bands separately and produced S/N of 4.2, 6.3 and 
5.5 for the WAPPs centered at 1410, 1510 and 1610 MHz. Although S/N = 9.4 is close to the detection 
threshold for pulsars in a blind search, it is much more significant in this case because of the reduced number 
of trials in this search relative to a blind search. 
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25-27 K o ff the plane of the G alaxy plus 6 K of Galactic emission in the specific direction of 
the pulsar lHaslam et al.lll982l ). Gain = 10.5 K Jy _1 and 2 polarizations, and an inefficiency 
factor of 12% due to the 3-level sampling of the WAPP correlators, we obtain for the first 
detection a flux density at 1.4 GHz of 5*1400 — 4.1 /iJy and for the second detection S140Q — 
3.1 /iJy. These values are consistent given the large relative uncertainties in the S/N estimates 
and the varying effect of radio frequency interference; at this DM, scintillation is not likely 
to cause a large variation in the flux density. 

The time- averaged flu x density is ~ 3.4 nJy. Using the NE2001 model for the electron 



distribution in the Galaxy (jCordes fc Lazioll2002l ). we ob tain from the pulsar's position and 
DM a distance of 3.2 kpc with a nominal error of 20% flCordes &: Laziol l2002l ). The time- 
averaged flux density thus corresponds to a pseudo-luminosity L M00 ~ 0.035 mJy kpc 2 . This 
is fainter than the least luminous young pulsar in the ATNF catalog (PSR J0205+6449, with 
a 1.4 GHz pseudo-luminosity of 0.5 mJy kpc 2 ). It is, however, more luminous than the radio 
pulsations discovered through a deep search of an other pulsar first di scovered by Fermi, PSR 
J1741— 2054 which has L1400 ~ 0.025 mJy kpc 2 (jCamilo et al.l 120091 ) . These two detections 
clearly demonstrate that some pulsars, as seen from the Earth, can have extremely low 
apparent radio luminosities; i.e., similarly deep observations of other 7-ray selected pulsars 
might detect additional very faint radio pulsars. We note that these low luminosities, which 
may well be the result of only a faint section of the radio beam crossing the Earth, are much 
lower than what has often been termed "radio quiet" in population s ynthesis models used t o 



estimate the ratio of "radio- loud" to "radio quiet" 7-ray pulsars (eg. iGonthier et al.l 12004 ) . 



185 4. Energy-Dependent Gamma-ray Pulse Profiles 

186 The pulse profile and spectral results reported in this paper use the survey data collected 

187 with the LAT from 2008 August 4 through 2009 September 18. We selected "diffuse" class 

188 photons (see j|2]) with energies E > 100 MeV and, to limit contamination from photons from 

189 Earth's limb, with zenith angle < 105°. 

190 To explore the dependence of the pulse profile on energy, we selected an energy- dependent 

191 region of interest (ROI) with radius 9 = 0.8 x E~ ' 75 degrees, but constrained not to be out- 

192 side the range [0.35°, 1.5°]. We chose the upper bound to minimize the contribution from 

193 nearby sources and Galactic diffuse emission. The lower bound was selected in order to in- 

194 elude more photons from the wings of the point spread function (PSF) where the extraction 

195 region is small enough to make the diffuse contribution negligible. Figure |2] shows folded 
we light curves of the pulsar in 32 constant-width bins for different energy bands. We use the 
197 centroid of the 1.4 GHz radio pulse profile to define phase 0.0. Two rotations are shown in 
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we each case. The top panel of the figure shows the folded light curve for photons with E > 0.1 

199 GeV. The 7-ray light curve shows two peaks, PI at phase 0.220 ± 0.002 which determines 

200 the offset with the radio peak, 5. The second peak in the 7-ray, P2, occurs at phase 0.580 

201 ± 0.003. The phase separation between the two peaks is A = 0.360 ± 0.004. The radio 

202 lead 5 and gamma peak separation A valu es are in good agreement wi th the correlation 



203 predicted for outer magnetosphe re models, (Romani fc Yadigaroglulll995l ) and observed for 



204 other young pulsars (Figure 3 of lAbdo et al. 



(|2009cJ)). 



205 Pulsed emission from the pulsar is clearly visible for energies E > 5 GeV with a chance 

206 probability of ~ 4 x 10~ 8 . Pulsed emission is detected for energies above 10 GeV with a 

207 confidence level of 99.8%. We have measured the integral and widths of the peaks as a 
20B function of energy and have found no evidence for significant evolution in shape or P1/P2 

209 ratio with energy. We note that the pulsar is at low Galactic latitude (b ~ —0.89° ) where 

210 the Galactic 7-ray diffuse emission is bright ( it has not been subtracted from the light curves 

211 shown.) 

212 Figure E] shows the observed LAT counts map of the region around PSR J1907+0602. 

213 We defined the "on" pulse as pulse phases 0.12 < <fi < 0.68 and the "off" pulse as its 

214 complement (0.0 < <f> < 0.12 and 0.68 < (f> < 1.0). We produced on-pulse (left panel) and 

215 off-pulse (right panel) images, scaling the off-pulse image by 1.27. The figure indicates the 

216 complexity of the region that must be treated in spectral fitting. Besides the pulsar there 

217 are multiple point sources, Galactic, and extragalactic diffuse contributions. 



5. Energy Spectrum 

The phase-averaged flux of the pulsar was obtained by performing a maximum likelihood 
spectral analysis using the Fermi LAT science tool gtlike. Starting from the same data set 
described in [gj we selected photons from an ROI of 10 degrees around the pulsar position. 
Sources from a preliminary versio n (based on 1 1 months of data) of the first Fermi LAT 



223 7— ray catalog f lAbdo et al.l l2009cf ) that are within 15 degree ROI around the pulsar were 

224 modeled in this analysis. Spectra of sources farther away than 5° from the pulsar were fixed 

225 at the cataloged values. Sources within 5° degrees of the pulsar were modeled with a simple 

226 power law. For each of the sources in the 5° degree region around the pulsar, we fixed the 

227 spectral index at the value in the catalog and fitted for the normalization. Two sources 

228 that are at a distance > 5° showed strong emission and were treated the same way as the 

229 sources within 5°. The Galactic diffuse emission (gll_iemv02) and the extragalactic diffuse 
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230 background (isotropic_iem_v02) were modeled as well. 

The assumed spectral model for the pulsar is an exponentially cut-off power law: dN/ dE - 
N Q (E / E )~ r exp(— E/E c ). The resulting spectrum gives the total emission for the pulsar 
assuming that the 7-ray emissi on is 100% pulsed. The unbinned gtlike fit, using P6_v3 



instrument response functions (lAtwood et al.l 120091 ). for the energy range E > 100 MeV 



gives a phase-averaged spectrum of the following form: 

dN 
~d~E 



% = (7.06 ± 0.45W. + aZUs.) x 10- u E- r e- E ^ 7 cm" 2 s" 1 MeV" 1 (1] 



231 where the photon index T = 1.76 ± 0.05 sta t. + (-o^l^sys. and the cutoff energy E c = 3.6 ± 

232 0.5 sta t. + (In ll)svs. GeV. 



The integrated energy flux from the pulsar in the energy range E > 100 MeV is F, 
(3.12 ± 0.15 stat + {tolDsys.) x 10~ 10 ergs cm~ 2 s _1 . This yields a 7-ray luminosity of L„ 



7 



235 AnfaF^d 2 = 3.8 x 10 35 /n^i.2 er g s s_1 above 100 MeV, where /q is an effective beaming factor 

236 and d 3 2 — d/(3.2)kpc. This corresponds to an efficiency of 77 = L 1 /E = 0.13F 7 d| 2 f° r 

237 conversion of spin-down power into 7— ray emission in this energy band. 

23B We set a 2a flux upper limit on 7-ray emission from the pulsar in the off-pulse part 

239 of F a ff < 8.31 x 10 -8 cm -2 s -1 . In addition to the 7-ray spectrum from the point-source 

240 pulsar PSR J1907+0602, we measured upper limits on 7-ray flux from the extended source 

241 HESS J1908+063 in the energy range 0.1-25 GeV. We performed binned likelihood analysis 

242 using the Fermi Science Tool gtlike. In this analysis we assumed an extended source 

243 with gaussian width of 0.3° and 7-ray spectral index of —2.1 at the location of the HESS 

244 source. The upper limits suggest that the spectrum of HESS J1908+063 has a low energy 

245 turnover between 20 GeV and 300 GeV. Figure @] shows the phase-averaged spectral energy 

246 distribution for PSR J1907+0602 (green circles). On the same figure we show data points 

247 from HESS for the TeV source HESS J1908+063 (blue circles) and the 2a upper limits from 

248 Fermi for emission from this TeV source. Figure [5] shows an off pulse residual map of the 

249 region around PSR J1907+0602. The timing position of the pulsar is marked by the green 

250 cross. The 5a contours from Milagro (outer) and HESS (inner) are overlaid. As can be seen 

251 from the residual map, there is no gamma-ray excess at the location of either the pulsar or 

252 the PWN. 



3 Descriptions of the models are available at 



http : / /fermi . gsfc . nasa. gov / 
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6. X-ray Counterpart 

A 23 ks ASCA GIS exposure of the EGRET source GeV J1907+0557 revealed an ~ 
8' x 15' region of possible extended hard emission surroundin g two point-like peaks lying ~ 15' 



to the southwest of PSR J1907+0602 ([Roberts et al.ll200ll ) and no other significant sources 



in the 44' ASCA FOV. A 10 ks Chandra ACIS-I image of the ASCA emission (ObsID 7049) 
showed it to be dominated by a single hard point source, CXOU J190718. 6+054858 with no 
compact nebular structure and just a hint of the several arcminute-scale emission seen by 
ASCA. CXOU J190718. 6+054858 seemed to turn off for ~ 2 ks during the Chandra exposure, 
suggesting that it may be a binary of some sort or else a variable extragalactic source. There 
is no obvious optical counterpart in the digital sky survey optical or 2MASS near infrared 
images, nor in a I band image taken with the 2.4m Hiltner telescope at MDM (Jules Halpern, 
private communication). This strongly suggests that it is not a nearby source. An absorbed 
power law fits the spectrum of this source well, with absorption n# = l-8l *g x 10 22 cm -2 
(90% confidence region), a photon spectral index Y = 0.9loi, and an average 2-10 keV flux 
of 4.4+5 8 x 10 _13 erg cm -2 s _1 (68% confidence region). The fit absorption is similar to the 
estimated total Galactic absorption from the H EASARC nH tool of 1.6 x 10 22 cm~ 2 based 



269 on the Dickey and Lockman (1990) HI survey (jDickey fc Lockmanl Il990l ). suggesting that 

270 an nu of ~ 2 x 10 22 cm~ 2 is a reasonably conservative estimate of interstellar absorption for 

271 sources deep in the plane along this line of sight. 

272 The timing position of LAT PSR J1907+0602 is in the central 20' of the ASCA GIS 

273 FOV (Figure [6]). There is no obvious emission in the ASCA image at the pulsar position. 



274 Using the methodology of iRoberts et al.l (120011 ). a 24 pixel radius extraction region (~ 6' 



275 and assuming an absorbed power law spectrum with nn = 2 x 10 22 cm -2 and T — 1.5, we 

276 place a 90% confidence upper limit on the 2-10 keV flux F x < 5 x 10 _14 erg cm -2 s" 1 . This 

277 suggests that for any reasonable absorption, the total unabsorbed X-ray flux from the pulsar 
27B plus any arcminute-scale nebula is less than 10 _13 erg cm" 2 s -1 . 

279 PSR J1907+0602 was well outside of the FOV of the first Chandra observation, and 

280 so we proposed for an observation centered on the pulsar. We obtained a 19 ks exposure 

281 with the ACIS-S detector (ObsID 11124). The time resolution of the ACIS-S detector on 

282 board Chandra does not allow for pulse studies. The only source within an arcminute 

283 of the timing position and the brightest source in the FOV of the S3 chip is shown in 

284 Figure It is well within errors of the timing position. Examination of the X-ray image 

285 in different energy bands showed virtually no detected flux below ~ IkeV and significant 

286 flux above 2.5 keV, suggesting a non-thermal emission mechanism for much of the flux. 

287 A comparison of the spatial distribution of counts between 0.75 keV and 2 keV to those 

288 between 2keV and 8keV shows some evidence for spatial extent beyond the point spread 
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function for the harder emission but not for the softer emission. This would be consistent 
with an interpretation as predominantly absorbed but thermal emission from a neutron star 
surface surrounded by non-thermal emissio n from a compact pulsar wind nebula, which is 
the typical situation for young pulsars (see iKaspi et al.l 120061 . and references therein). We 
plot the Chandra 0.75-2keV, 2-8keV, and 0.75-8 keV images with an ellipse showing the 
timing position uncertainty, and a circle with a radius of 0.8". From a modeled PSF, we 
estimate 80% of the counts should be contained within this circle. While this seems to be 
the case for the soft image, only roughly half the counts in the harder image are contained 
within that radius. With only ~ 12 source counts in the 0.75-2 keV image within 6" and 
~ 30 source counts in the 2-8 keV image, quantitative statements about source size and 
spectrum are difficult to make. We obtain a best fit position for the nominal point source of 
R.A. = 19 h 07 m 54 f !76, decl. = +06°02'14.6" and estimated error of 0.7" (an additional 0.1" 
centroid fitting uncertainty added to the nominal Chandra 0.6" uncertainty). Using a 6" 
radius extraction region and an annulus between 6" and 24" for background, we extracted 
source and background spectra and fit them within XSPEC (Figure [8]). A simple power law 
plus absorption model fit the data well in the energy range 2-10 KeV, with best fit values 
n H = 1.3 x 10 22 cm" 2 and T = 1.6, with a total flux of 2.3t?;^ x lO^ergcm-V 1 . The 
low count rates and covariance between the absorption and photon index meant the spectral 
parameters could not be simultaneously meaningfully constr ained. Fixing the spectral index 
T = 1.6, a typical value for compact pulsar wind nebulae (IKaspi et al.l 120061 ). we obtain a 
90% confidence region for the absorption of 0.7 — 2.5 x 10 22 cm -2 , consistent with a source a 
few kilo parsecs or more away and with CXOU J190718. 6+054858 discussed above. We note 
that with such an absorption a significant thermal component in the below 2 keV emission 
is neither required nor ruled out by the spectral fitting. 



Discussion 



314 
315 
316 
317 
318 
319 
320 
321 
322 
323 



The dispersion measure from the radio detection suggests a distance of 3.2 kpc, with 
a nominal error of 20%. However, there are many outliers to the DM error distribution, 
although the l argest fractional errors tend to be from pulsars at high Galactic latitudes or 
very low DMs (beller et al. l l2009[ Ichatteriee et~alll2009h . For PSR J1907+0602, at a latitude 
b = —0.9° with a moderate DM, the distance estimate is likely to be reasonable. Since the 
apparent 7-ray pulsed efficiency in the Fer mi pass-band is well above the median for other 
gamma-ray pulsars in lAbdo et al.l (j2009dl ) (13% compared to 7.5%), it is worth checking 
secondary distance indicators to see if the DM measure could be a significant overestimate 
of the true distance. We can use the X-ray observations of PSR J1907+0602 to do this. 
Several authors have noted a correlation between the X-ray luminosity of young pulsars and 
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324 their spin-down power (eg. Saito 1998, Possenti et al. 2002, Li, Lu and Li 2007). Most 

325 of these have the problem of using X-ray fluxes derived from the literature using a variety 

326 of instruments with no uniform way of choosing spectral extraction regions. This can be 

327 especially problematic with Chandra data, since faint, arcminute scale emission can easily be 
32B overlooked. We compare our ASCA GIS upper limits to Figure 1 of Cheng, Taam and Wang 

329 (2004) who used only ASCA GIS data to derive their X-ray luminosity relationships. We see 

330 that the typical X-ray luminosity in the ASCA band for a pulsar with E = 2 x 10 36 erg s _1 

331 is L x ~ 10 33 — 10 34 ergs _1 with all of the pulsars used in their analysis with E > 10 36 ergs _1 

332 having L x > 10 32 ergs _1 . From these values and the ASCA upper limit, we derive a lower 

333 limit for the distance to LAT PSR J1907+0602 of ~ 3 kpc. From Figure 2 of Li, Lu and 

334 Li (2007), who used XMM-Newton and Chandra derived values, we see we can expect the 

335 luminosity to be between ~ 10 3L5 — 10 34,5 erg s _1 . From our detection with Chandra, we 

336 again estimate a lower distance limit of ~ 3 kpc. The "best guess" estimate from their 

337 relationship would result in a distance of ~ 13kpc. We note that if we assume the pulsed 
33B emission to be apparently isotr opic (i.e. fa = 1 as s imple outer gap models suggest should 



339 approximately be the case, see IWatters et al.l (120091 )). a distance of 9 kpc would result in 

340 100% 7-ray efficiency 

341 The derived timing position of PSR J1907+0602 is well inside the extended HESS source, 

342 although ~ 14' southwest of the centroid. The TeV source is therefore plausibly the wind 

343 nebula of PSR J1907+0602. The physical size of this nebula is then ^ 40 pc, and the 

344 integrated luminosity above 1 TeV is ^ 40% that of the Crab, and in the MILAGRO band 

345 (~ 20 TeV) at least twice that of the Crab. There is a hint of some spatial dependence of 

346 the TeV spectrum in the HESS data, with the harder emissi on (> 2.5 TeV) peaking nearer 



347 the pulsar than the softer emission ( lAharonian et al.ll2009l ). If confirmed, this would be 



34B consistent with the hardening of the TeV emis sion observed towards P SR B1823— 13, thought 



349 to be the pulsar powering HESS J1825— 137 ( lAharonian et al.ll2006l ). This latter pulsar has 

350 a spin period, characteristic age, and spin-down energy similar to PSR J1907+0602, and 

351 is also located near the edge of its corresponding TeV nebula. We also note that HESS 

352 J1825— 137 subtends ~ 1° on the sky and has a flux level above 1 TeV of around 20% of the 

353 Crab. While the overall spectrum of HESS J1825— 137 is somewhat softer than the spectrum 

354 of HESS J1908+063, near the pulsar its spectrum is similarly hard. At a distance of ~ 4 kpc, 

355 HESS J1825— 137 has a luminosity similar to the Crab TeV nebula, but with a much larger 

356 physical size of ~ 70pc. Given the distance implied above and a flux above 1 TeV ~ 17% of 

357 the Crab, HESS J1908+063 is similar in size and luminosity to HESS J1825— 137. 

35B At 20 TeV, HESS J1908+063 has a flux ~80% of the Crab, and so at a distance ^ 1.5 

359 times that of the Crab, is much more luminous at the highest energies. This is because there 

360 is no sign of a high-energy cutoff or break, as is seen in many other TeV nebulae. Aharonian 
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361 et al. (2009) place a lower limit of 19.1 TeV on any exponential cutoff to the spectrum. 

362 This implies that either th e spec trum is uncooled due to a very low nebular magnetic field 

363 ( 3/iG, see, eg. Ide Jagerl (120081 )). an age much less than the characteristic age of 19.5 kyr, 

364 or else there is a synchrotron cooling break below the HESS band. 

365 Our upper limits above a few GeV (Figure H|) requires there to be a low energy turnover 

366 between 20 GeV and 300 GeV. Given the nominal PWN spectrum, we constrain the overall 

367 PWN flux to be < 25% of that of the pulsar. If only the HESS band is considered, and 

368 assuming the DM distance, the TeV luminosity L PWN = 5 — 8%E. However, since the TeV 

369 emission is generally thought to come from a relic population of electrons the luminosity is 

370 likely a functio n of the spin-d own history of the pulsar rather than the current spin-down 



37i luminosity (eg. Ide Jagerl 120081 ). These numbers support consistency of the association of the 



372 TeV source with the pulsar, in the weak sense of not being discrepant with other similar 

373 systems. 



7.1. On the possible association with SNR G40.5— 0.5 



The bulk of HESS J1908+063 is between PSR J1907+0602 and the young radio SNR 
376 G40.5— 0.5, sug gesting a possible as sociation. The distance estimate (~ 3.4 kpc I Yang et al. 



20061 ) and age (IDownes et all Il980l ) estimates of SNR G40.5— 0.5 are also consistent with 



378 th ose of PSR J1907+060 2. If we use the usually assumed location for SNR G40.5-0.5 given 

379 by lLangston et al.l ( 120001 ) (RA=19 h 07 m lr!9, Dec=6°35'15"), we get an angular separation of 

380 ~ 35' between the timing position for the pulsar and the SNR. However, this position for the 

381 SNR is from single dish observations that were offset towards one bright side o f the nominal 

382 shell. We use the VLA Galactic Plane Survey 1420 MHz image ( jStil et al.ll2006l ) of this region 

383 to estimate the SNR center to be RA=19 h 07 m 08!6, Dec=6°29'53" (Figure EJ) which, for an 

384 assumed distance of 3.2 kpc, would give a separation of ~28 pc. Given the characteristic age 

385 of 19.5 kyr years, this would require an average transverse velocity of ~ 1400 km/s. While 

386 velocities about this high are seen in some cases (eg. PSR B 1508+55 has a transverse velocity 

387 of ~ 1100 km/s, Chatterjee et al. 2005 ), it is several times the average pulsar velocity and 

388 many times higher than the local sound speed. We note that pulsars with a braking index 

389 significantly less than n = 3 assumed in the d erivation of the ch aracteristic age could have 

390 ages as much as a factor of two greater (see eg. iKaspi et al.l 1200 ll ). and thus a space velocity 

391 around half the above value may be all that is required. But with any reasonable assumption 

392 of birthplace, distance, and age, if the pulsar was born in SNR G40.5— 0.5, any associated 

393 X-ray or radio PWN should show a bow-shock and trail morphology, with the trail likely 

394 pointing back towards the SNR center. Unfortunately, the compactness and low number of 
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counts in our Chandra image precludes any definite statement about the PWN morphology. 
One arrives at a different, and lower, minimum velocity if one assumes the pulsar was born 
at the center of the TeV PWN and moved to its present position, but the resulting velocity 
would still require a bow shock. 

One can also get a pulsar offset towards the edge of a relic PWN if there is a significant 
density gradient in the surrounding ISM. A gradient will cause the supernova blast wave to 
propagate asymmetrically. Where the density is higher, the reverse shock propagating back 
to the explosion center will a lso be asymmetric. This will tend to push the PWN away from 
the region of higher density (IBlondin et al.l l200ll ; iFerreira &: de Jagerl 120081 ) . This has been 



invoked to explain the offsets in the Vela X and HESS J1825— 137 nebulae as well as several 
others. Infrared and radio imaging of the region shows that HESS J1908+063 borders on 
a sh ell of ma t erial surrounding the S74 HII region, also known as the Lynds Bright Nebula 
352. iRusseill (120031 ) gives a kinematic distance of 3.0 ± 0.3 kpc for this star forming region, 
compatible with the pulsar distance. In this scenario, the pulsar would not have to be highly 
supersonic to be at the edge of a relic nebula, and would not have to be traveling away from 
the center of the TeV emission. 

A third, hybrid possibility is that SNR G40.5— 0.5 is only a bright segment of a much 
larger remnant, whose emission from the side near the pulsar is confused with that from the 
molecular cloud. The asymmetry would be explained by the difference in propagation speed 
in the lower density ISM away from the molecular cloud. 

Our current Chandra data are insufficient to distinguish between the above scenarios. 
However, there is also the possibility of a compact co metary radio nebula, such as is seen 



aroun d PSR B1853+01 in SNR W44 (IFrail et al.lll996f ) and PSR B0906-49 (IGaensler et al. 



19981 ). In addition, sensitive long wavelength radio imaging could reveal any larger, faint 
SNR shells. Imaging with the EVLA and LOFAR of this region is therefore highly desirable. 
The connection between the pulsar and the TeV nebula could be further strengthened by a 
confirmation of the spatio-spectral dependence of the nebula where the spectrum hardens 
nearer to the pulsar. 
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Table 1: Measured and Derived timing parameters of PSR J1907+0602 



Fit and data-set 

Pulsar name J1907+0602 

MJD range 54647-55074 

Number of TOAs 23 

Rms timing residual (/is) 404 

Measured Quantities 

Right ascension, a 19:07:54.71(14) 

Declination, 5 +06:02:16.1(23) 

Pulse frequency, v (s" 1 ) 9.3780713067(19) 

First derivative of pulse frequency, i> (s -2 ) . . —7.6382(4) xl0 _1 

Second derivative of pulse frequency, v (s~ 3 ) 2.5(6) xl0~ 22 

Epoch of frequency determination (MJD) . . 54800 

Dispersion measure, DM (cm _3 pc) 82.1(11) 

Derived Quantities 

Characteristic age (kyr) 19.5 

Surface magnetic field strength (G) 3.1 x 10 12 

E (erg s" 1 ) 2.8 x 10 36 

Assumptions 

Time units TDB 

Solar system ephemeris model DE405 



Note. — The numbers in parentheses are the errors in the last digit of the fitted parameters. The errors 
are statistical only, except for the position error, as described in j Ql The derive d parameters of E, B, and r c 
are essentially unchanged with respect to those reported by ([Abdo et al]|2009af) , but the position has moved 
by 1.2 arcmin. 
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Fig. l.- 
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Post-fit timing residuals for PSR J1907+0602. The reduced chi-square of the fit is 
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Fig. 2. — Folded light curves of PSR J1907+0602 in 32 constant-width bins for different 
energy bands and shown over two pulse periods with the 1.4 GHz radio pulse profile plotted 
in the bottom panel. The top panel of the figure shows the folded light curve for photons 
with E > 0.1 GeV. The other panels show the pulse profiles in exclusive energy ranges: 
E > 3.0 GeV (with E > 5.0 GeV in black) in the second panel from the top; 1.0 to 3.0 GeV 
in the next panel; 0.3 to 1.0 GeV in the fourth panel; and 0.1 to 0.3 GeV in the fifth panel. 
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Fig. 3. — The observed Fermi-LAT counts map of the region around PSR J1907+0602. Left: 
"on" pulse image, right: "off" pulse image. The open cross-hair marks the location of the 
pulsar. Color scale shows the counts per pixel. 
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Fig. 4. — Phase-averaged spectral energy distribution for PSR J1907+0602 (green circles). 
Blue circles are data from HESS for HESS J1908+063 TeV source. 2 a upper limits from 
Fermi for emission from this TeV source are shown in blue. The black line shows the spectral 
model for the pulsar (equation CD). The upper limits suggest that the spectrum of HESS 
J1908+063 has a low energy turnover between 20 GeV and 300 GeV. 
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Fig. 5. — Residual map of the region around PSR J1907+0602 in the off-pulse. The timing 
position of the pulsar is marked by the cross. The 5 a contours from Milagro (outer) and 
HESS (inner) are overlaid. 
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Fig. 6.— ASCA GIS 2-10 keV image of the region around PSR J1907+0602. The green 
contours are the 4-7 a significance contours from HESS. 
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Fig. 7. — Chandra ACIS images of PSR J1907+0602. The blue ellipse shows the uncertainty 
in the timing position. The green circle of radius 0.8" is twice the FWHM of the 5keV PSF 
at this position, and should contain roughly 80% of the counts. The image at 0.75-2 keV 
(Left), 2-8 keV (Center) and 0.75-8 keV (right) is shown. Color scale shows the counts per 
pixel. 
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Fig. 8.- 



Chandra X-ray spectrum of PSR J1907+0602. 
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Fig. 9. — VGPS 1420 MHz image of region in Galactic coordinates showing relationship 
between SNR G40.5— 0.5, HESS J1907+063 (blue contours representing the 4,5,6 and la 
significance levels), the star forming region S74, and PSR J1907+0602. 



